A chiral lagrangian with Broken Scale: 
testing the restoration of symmetries in astrophysics and in the laboratory 
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We study matter at high density and temperature using a chiral lagrangian in which the breaking 
of scale invariance is regulated by the value of a scalar field, called dilaton [l], 0,0,0] . We provide a 
phase diagram describing the restoration of chiral and scale symmetries. We show that chiral sym- 
metry is restored at large temperatures, but at low temperatures it remains broken at all densities. 
We also show that scale invariance is more easily restored at low rather than large baryon densities. 
The masses of vector mesons scale with the value of the dilaton and their values initially slightly 
decrease with the density but then they increase again for densities larger than ~ 3po- The pion 
mass increases continuously with the density and at po and T=0 its value is ~ 30 MeV larger than in 
the vacuum. We show that the model is compatible with the bounds stemming from astrophysics, as 
e.g. the one associated with the maximum mass of a neutron star. The most striking feature of the 
model is a very significant softening at large densities, which manifests also as a strong reduction of 
the adiabatic index. While the softening has probably no consequence for Supernova explosion via 
the direct mechanism, it could modify the signal in gravitational waves associated with the merging 
of two neutron stars. 



PACS numbers: 21.65.Mn, 12.39.Fe, 21.65.Cd, ll.10.Wx 



o 

=5 



(N 
> 
00 

oo 



in 
o 

oo 
O 



I. INTRODUCTION 

A very hot problem in hadronic and nuclear physics 
is to investigate the Equation of State (EOS) of mat- 
ter at large density and/or temperature. Several exper- 
iments have been proposed and will likely be performed 
in the next years. One experiment has been planned at 
facility FAIR at GSI [5| , other experiments could be per- 
formed at RHIC (Brookhaven) and at the Nuclatron in 
Dubna. In all these experiments Heavy Ion Collisions 
(HICs) will take place at energies of the order of a few 
ten A GeV, which correspond to the regime at which 
maximum baryon densities can be reached. Another im- 
portant source of information concerning matter at high 
densities is the study of the astrophysics of compact stars, 
where the temperatures are moderate (typically lower 
than a few ten MeV), but the baryon density can reach 
several times nuclear matter saturation density po. It 
is therefore extremely important to investigate the EOS 
by using various theoretical techniques, providing hints 
about the main features of matter at those extreme con- 
ditions. In a sense one is providing a map of the "new" 
regions which experiments will likely explore. 

In this work we use an effective model for nuclear 
physics, as a tool for exploring theoretically the unknown 
regions. The main feature of this model is that it is based 
on a chiral invariant lagrangian. It is well known that it 
is not trivial to introduce chiral symmetry in an effec- 
tive nuclear lagrangian, so that saturation properties of 
nuclei are well reproduced together with the possibility 
to restore chiral symmetry at large density and/or tem- 
perature. For instance, the attempt of using the linear 
sigma model to describe nuclear dynamics fails due to 
the impossibility of reproducing the basic properties of 
nuclei Q. The coefficients of the self-couplings of the 



scalar field are dictated by the form of the potential and 
chiral symmetry restoration takes place already at about 
Pq. More sophisticated approaches have been proposed 
in the literature. A few works are based on a SU(2) chi- 
ral symmetry scheme, either adopting the parity-doublet 
formulation 0, d, @ or introducing the scale invariance 
[Tol[ilT|. O ther works extend the symmetry to the strange 
sector [12l . H3I [bH . [lol [r| . Here we will explore the model 
introduced by the nuclear physics group of the University 
of Minnesota [E 0j B 0| • I n that model chiral fields are 
present together with a dilaton field which reproduces 
the breaking of scale symmetry in QCD, but at a mean 
field level. In their first paper [l[ they have shown that 
it is indeed possible to reproduce the main features of 
closed shell nuclei; in a second paper [2| the phenomenol- 
ogy of chiral fields, including their interaction with the 
nuclei, is well described; in a third paper [3| they discuss 
mesons at finite temperature showing in particular that 
chiral symmetry is restored at large temperature; finally 
in Ref. [J] they discuss symmetric nuclear matter at finite 
density and finite temperature. 

In our paper we extend the previous analysis providing 
rather detailed previsions concerning various quantities 
which can be measured in the laboratory, like e.g. the 
compressibility of matter and the behavior of the masses 
of meson fields at large densities. We will study also 
the dependence of those quantities on the isospin frac- 
tion. Moreover, we will test this chiral model at large 
baryon and isospin densities which are relevant for study- 
ing various astrophysical problems, as the structure and 
formation of compact stars. In this way we complete 
the analysis provided in Ref. .17| where we compared the 
compressibility computed in various nuclear and hadronic 
models for the EOS of matter. 

A few caveats concerning the validity of our results 
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at very large densities and temperatures are in order: 
clearly hyperonic and/or quark degrees of freedom should 
appear in the high energy density regime. On the 
other hand before investigating even more complicated 
lagrangians it is worth analyzing in detail the outcome 
of a purely hadronic model. It is also important to re- 
mark that the restoration of scale invariance was ana- 
lyzed only in Refs. d [j] obtaining a unrealistically large 
critical temperature. In our paper we show that the in- 
clusion of thermal fluctuations of the vector meson fields 
reduces the critical temperature for scale restoration to 
a more realistic value and provides a rather new and in- 
teresting phase diagram. Hyperons are not included in 
the present work in order to keep the discussion of chi- 
ral symmetry and scale invariance restoration as simple 
as possible. A more realistic model should incorporate 
also h ype rons, following e.g. the schemes proposed in 
Refs. |ll [H, Q [H Notice finally that the la- 

grangian we are investigating mimics the behavior of the 
gluon condensate through the dynamics of the dilaton 
field. Since the value of the gluon condensate can de- 
termine if the system settles in the hadronic or in the 
quark-gluon phase [l8j , this model is particularly suit- 
able to be extended to the quark-gluon sector in order to 
investigate deconfinement, what will be done in a future 
work. 

The structure of the paper is the following. In Sec. |TT] 
we will describe the model we are using, in Sec. Mil we 
discuss the EOS as it can be tested in lab experiments, 
in Sec. IIVI we show the implications for astrophysics. Fi- 
nally, in Sec. [V] we will present our conclusions. 



II. THE CHIRAL-DILATON MODEL 



chosen so that it reproduces (at mean field level) the 
divergence of the scale current which in QCD is given 
by eq. (pj. The potential of the dilaton field is therefore 
determined by the equation: 

^ = 4^)-^=4e TOC (£)\ (2) 

where the parameter e vac represents the vacuum energy. 

The previous approach holds in the case of pure gauge 
QCD where, at first loop, the violation of the scale in- 
variance reads: 

where N c is the number of colors. To take into account 
massless quarks a generalization was proposed in [2l| , so 
that also chiral fields contribute to the trace anomaly. 
In this way the single scalar field of eq. ((2]) is replaced 
by a set of scalar fields {a, n ,</)}. The relative weight 
of the gluons and of the quarks to the violation of scale 
invariance can be read from the QCD /3-function, which 
at first loop is: 

where n/ is the number of massless flavors. 

B. The lagrangian 

The lagrangian of the chiral-dilaton model (CDM) 
reads: 



A. Scale invariance 

The aim of many effective lagrangians for nuclear 
physics is to mimic the behavior of QCD in the strong 
regime, where chiral and scale symmetries are broken. 
While the breaking of chiral symmetry is connected with 
the generation of fermionic masses, the breaking of scale 
invariance is due to quantum corrections (the so-called 
scale or trace anomaly) and is likely associated with the 
formation of a gluon condensate. More explicitly, the 
violation of scale invariance in QCD corresponds to a 
non-conservation of the dilatation current, which reads: 



(1) 



where Jqcd * s the dilatation current in QCD, 9^ is 
the "improved" energy- momentum tensor, F£ v (x) is the 
gluon field strength tensor and /3(g s ) is the QCD beta 
function. 

In the approach of Schechter, Migdal and Shifman [l^, 
|20| a scalar field representing the gluon condensate is 
introduced and its dynamics is regulated by a potential 
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where the potential governing the chiral and the dilaton 
field reads: 
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Here a and n are the chiral fields, <j> the dilaton field, 
uj^ the vector meson field and b M the vector-isovector p- 
meson field, introduced here in order to study asymmet- 
ric nuclear matter. The field strength tensors are defined 
in the usual way = d^Wv—dvW^ B Mt/ = d^h^-d^h^. 
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The vacuum mean field value of the scalar fields is: 
4> — 4>o, a — ct an d 7f = 0. Moreover, C = ( t>o/ cr Oi B 
and 5 are constants and is a term which breaks explic- 
itly the chiral invariance of the lagrangian. 

It is easy to check that the first two terms of the po- 
tential V and the chiral symmetry breaking terms con- 
tribute to the trace anomaly, so that eq. $2$ is satisfied 
with e vac = —jB<fio' 1 (l — 5) — e' 1 , while the third term fixes 
the vacuum conditions. Moreover, since e vac takes contri- 
butions both from gluons and quarks, the choice 8 = 33 
is suggested by the form of the QCD beta-function at 
first loop, Eq. Q. 

Notice that, as in the case of the well known "Mexican 
hat" potential, the minimum of the potential V is degen- 
erate if d x — and the ground state satisfies the "chiral 
circle" equation (ct 2 +7r 2 )/er 2 = 4 >2 /4'o 2 - Notice there- 
fore that when the mean value of the dilaton field drops 
to zero also the radius of the chiral circle drops to zero. 
In other words, the restoration of scale invariance implies 
the vanishing of all the condensates, including the chiral 
one. 

It is important to remark that scale invariant masses 
for the vector mesons can be generated by coupling the 
vector fields either to the dilaton or to the chiral fields. 
In the initial versions of the model [U, [22j the coupling 
\ \_G uja {a 2 + 7r 2 )] ui^lo^ was adopted, which leads to the 
scenario proposed by Brown and Rho where the vector 
meson masses scale as the mass of the nucleon (23j . As 



TABLE I: Values of the parameters. 



it will be clarified in the following, that choice leads to 
an anomalous behavior of the sigma mean field which in- 
creases at very large densities. In Ref. [l[ a more gen- 
eral combination \ [G^^tp 2 + G LUa (a 2 + 7r 2 )] lu^uj 11 was 
tested, but that choice does not provide a good phe- 
nomenology for nuclei. In the lagrangian of Eq. ^ the 
uj and p masses are generated by coupling the vector 
mesons only to the dilaton field so that in the vacuum 
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and 



1/2 

= GA (pa. This choice, success- 
fully adopted in Refs. @, H, does not lead to the sce- 
nario proposed by Brown and Rho [23[ . 

The values of the parameters which we used for cal- 
culations are listed in Table HI They were determined in 
Ref. [IJ by fitting the properties of nuclear matter and 
finite nuclei. In Ref. t 4] two parameter sets are provided, 
one with G4 — and the other with G4 ^ 0. Since a 
better phenomenology was obtained with the non-zero 
value of G4, we use that parameter set in our calcula- 
tions. The value of the coupling constant of the p-meson 
to the nucleon g p has been chosen in order to obtain for 



the symmetry energy at po the value E s 



35 MeV. 



The p-meson did not appear in Ref. [4] where only isospin 
symmetric nuclear matter has been studied. 
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1.41 


cr (MeV) 


102 


e'x 1/4 (MeV) 


119 



masses of the meson fields it is mandatory to take into 
account the thermal fluctuations of those fields. In the 
present work we consider the thermal fluctuations of the 
chiral fields and of the vector mesons. The fluctuations of 
the dilaton field are not taken into account in the present 
paper, due to the large mass of the glueball (about 1.6 
GeV). In a future analysis we will consider also the ther- 
mal fluctuations of the dilaton field, which become im- 
portant when studying in detail the temperature-density 
region where scale invariance is restored. 

To take into account thermal fluctuations one has to 
resort to a technique (developed in Refs.[3|, |4[) which al- 
lows to deal with non-polynomial functions of the fields. 
The technique is the following: 

a) the chiral invariant combination of the chiral fields 
is expanded in mean field value plus fluctuations, as: 
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where a is the sigma mean field value and Act is its fluc- 
tuation. The pion field has a vanishing mean field value 
and therefore 7r 2 is also the fluctuation squared; 

b) in order to simplify the evaluation of the thermal 
averages of the fluctuations, a crucial assumption is made 



(Act 2 



(8) 



This approximation is based on the observation that 
at low temperature the fluctuations are small (and 
therefore the error in assuming the validity of eq. ijHJl is 
also small), while at large temperature chiral symmetry 
is (at least partially) restored and eq. flSJ is a sensible 
approximation. We end up therefore with a unique 
quantity describing all the thermal fluctuations of the 
chiral fields, (ip 2 ), where xp 2 = (Act 2 



^ 2 )K 2 ; 



C. Thermal fluctuations and their resummation 

The main purpose of our work is to study the EOS 
of nuclear matter at finite temperature. When the tem- 
perature is large enough to become comparable to the 



c) a generic function of the chiral invariant expression 
is expanded as a series of the field fluctuations and the 
fluctuations in xp 2 are taken to all orders [651 ] : 

d) to evaluate the numerical value of {xjj 2 ) one needs 
the thermal averages of the fluctuations squared of the 
chiral fields, which are computed as integrals of the sta- 
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e) finally, the series in (ij> ) is resummed. 

In a more recent paper [24| a new technique to eval- 
uate thermal fluctuations has been developed, in which 
the fluctuations of a and tt are always considered as in- 
dependent quantities. In particular, a generic function of 
a and 7r is first expanded as a series of the chiral fields 
fluctuations and these fluctuations are taken to all or- 
ders. The series in (Act 2 ) and (tv 2 ) is then resummed. 
In the case of isospin symmetric matter the fluctuations 
of the three components of the pionic field are equal and 
the following result is obtained: 

/OO OO 
dzP a (z) / dyy 2 P n (y)f(a+ Z) y 2 ) 
-oo JO 

(11) 

where 

P„(*) = (2,r<AOr"»oxp (-27^3)) 

^(»)-v / !((^t) s,2 °' , (-#>) (12) 

are the gaussian weighting functions. This procedure can 
be generalized to the case of asymmetric matter, where 
the fluctuations of the pionic field are independent. In 
that case instead of the bidimensional integral of Eq. (fTTj) 
a four dimensional integral is obtained. Notice that the 
procedure developed in Ref. [HJ reduces to the technique 
discussed in Refs. [!, 0] if the thermal fluctuations of the 
chiral fields are assumed to be all equal. 

In the present work we have mainly adopted the 
technique of Refs.[l, 0| because our main aim is to 



extend the calculations towards asymmetric nuclear 
matter and astrophysics and we like to be able to 
directly compare our results with the ones obtained in 
previous calculations. Most of the results discussed in 
the following are rather independent from the technique 
adopted: e.g. the value of the sigma mean field is the 
same in both cases. Nevertheless it is necessary to use 
the technique developed in Ref.(24| when computing the 
isospin splitting of quantities strongly dependent on the 
thermal fluctuations, in particular the isospin splitting 
of the mass of the pion. 



D. Thermodynamics 

The value of the fields and of the various thermody- 
namical quantities at finite temperature and density can 
be obtained in two ways: i) by solving the field equations 
(obtained by minimizing the action) and substituting the 
solutions into the thermodynamical functions or, ii) by 
minimizing directly the relevant thermodynamical quan- 
tities, such as the free energy (if working at given density 
and temperature), or the thermodynamical potential (if 
working at given chemical potential and temperature) or 
the e nerg y density (if working at given entropy and den- 
sity) [251 ] . These two approaches are in principle equiv- 
alent, although the approximations discussed in the pre- 
vious subsection can introduce some (small) discrepancy 
4] . In the present paper we have decided to minimize di- 
rectly the thermodynamical quantities because in future 
works we aim at introducing quark and gluon degrees of 
freedom. Since in Ref. [H[ the pure gluonic sector was de- 
scribed in terms of an effective thermodynamical poten- 
tial (without introducing a model lagrangian) , the second 
approach is the appropriate one if quarks and gluons are 
introduced following the scheme proposed in [181 ]. 

The thermodynamical potential, within the mean field 
approximation and taking into account the thermal fluc- 
tuations, reads: 



A;\ 1 ^r, 2 J 2 , 2 1 rr , 2 J 2 h 2 r< 4 , A 1 ^* 2 /A^ 2 \ y^* 2 i^. 2 \ 
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where the quantities into brackets are thermal averages. 



We omit to show explicitly the complicate expression for 
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(V) , which can be obtained using the procedure described 
in the previous subsection and is given in Eq. (21) of 
Ref. [|| . The subscripts n and p stay for proton and neu- 
tron, respectively, the index a runs on the three isospin 
components of the pion and of the p-field and, finally, 
X = 0/<f>o- 



The energies which appear into the thermodynamical 
integrals read e* = J m* 2 + k 2 where m* is the effective 
mass of the i£/i-field, which can be obtained by comput- 
ing the thermal average of the second derivative of the 
lagrangian with respect to the field: 



M* N 


= gW* 2 
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2 2 

= rn^x , 
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m p 


2 2 

= m p X , 
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Va 2 + n 2 (a 2 +7T 2 ) 3 / 2 



V<7 2 + TT 2 (a 2 +7T 2 )3/ 2 , 



(14) 
(15) 
(16) 

(17) 
(18) 



where ps P and ps n are the scalar densities of the protons 
and of the neutrons, respectively. As already remarked 
above, we have omitted the thermal average on the dila- 
ton field because its thermal fluctuations are suppressed 
due to the large mass of this held. Notice that the masses 
of the chiral fields depend on their thermal fluctuations, 
which in turn depend on the masses. Due to these com- 
plicate dependencies, the minimization of the relevant 
thermodynamical quantities has to be done iteratively. 

Finally, the chemical potentials entering the thermo- 
dynamical functions are the effective ones. In the case of 
the nucleons, the effective chemical potentials are related 
to the standard ones as follows: 

Mp = Mp ~ 9w^o — \gpbo 

/< = Mn - g^Q + \g P h) ■ (19) 

Notice that, when asymmetric nuclear matter is consid- 
ered, the charged isovector mesons acquire a non zero 
chemical potential, due to the conservation of the isospin 
charge. The chemical potentials arc determined by the 
equations: 

Mtt+ = Mp+ = Mp ^ M« (20) 

and the effective ones, which enter e.g. Eq. ([13]) . are 
given by the relation: 

Pl+ = Mp+ = Mp - M,* = Mp - Mn - g P bo ■ (21) 

E. Goldstone theorem at finite temperature 

When discussing a field theory at finite temperature, 
one has to choose an approximation scheme which do not 



violate the symmetries of the lagrangian. In particular, 
adopting a chiral lagrangian and working in the chiral 
limit, at T=0 the pion is massless due to the Goldstone 
theorem. At finite temperature though, it is not auto- 
matic that the Goldstone theorem is satisfied, since a 
wrong truncation can spoil the theorem. There are sev- 
eral contributions to the various observables which one 
can in principle take into account: 

a) mean field contributions, neglecting all fluctuations; 

b) fluctuations in 7r 2 and u 2 , particularly at finite tem- 
perature; 

c) fermionic contributions associated with the so-called 
particle- hole excitations; 

d) further corrections, as the exchange diagrams dis- 
cussed e.g. in Ref. [13] at the end of Sec. IVA. 

A important guideline in deciding which truncation 
scheme is appropriate (in other terms, a way for decid- 
ing that including a further correction is not necessarily 
a better approximation) is to check what happens to the 
Goldstone theorem at finite temperature. For instance, 
in [3[ (which describes the same approximation scheme 
adopted here) the Goldstone theorem is satisfied also at 
finite temperature by taking into account contributions 
a) and b) and neglecting the others (see Fig. 3 of Ref. [|[ 
and Fig. [3] of the present paper). It is well known that, 
when studying chiral models, another truncation exists 
which allows the Goldstone theorem to be satisfied at fi- 
nite temperature and it corresponds to the inclusion of 
contributions a) and c) (see e.g. 25]). In the chiral limit 
one can check that indeed the contributions a) and c) 
to the mass mutually cancel. On the other hand, if one 
includes a) and b) and c), the Goldstone theorem is not 
satisfied at finite temperature, as already discussed in 
Ref. [13. 

In conclusion, one has to resort either to a) + b) (as we 



are doing here) or to a) + c) (as discussed e.g. in [25|). 
If one wants to add c) to a) and b), other contributions, 
as the ones in d) has also to be considered in order to 
satisfy Goldstone theorem. 

We will come back to this problem in Sec lIVDl when 
discussing the possibility that the pions form a Bosc- 
Einstcin condensate. 



III. MATTER AT HIGH BARYON DENSITY 
AND TEMPERATURE 

In this section we will focus on nuclear matter cal- 
culations at rather high temperature and relatively low 
isospin asymmetry. These calculations are aimed at pro- 
viding previsions on the EOS of matter in regions of the 
energy-density vs. baryon-density plane which can be 
tested with HIC experiments. Those experiments are 
particularly interesting at energies of the order of a few 
ten A GeV, which correspond to the regime at which 
maximum baryon densities can be reached. In this Sec- 
tion we first discuss the restoration of chiral symmetry 
and of scale invariance in various region of the energy- 
density vs. baryon density plane, we then analyze various 
phenomenological implications of the partial restoration 
of these two symmetries. 

A. Phase diagrams 

The advantage of using a sophisticated model which 
includes both the chiral dynamics and the gluon conden- 
sate dynamics, is that one can make previsions on the 
restoration of both the chiral and the scale symmetry. 
These previsions are different if one considers the case 
of exact chiral symmetry or the case where a symme- 
try breaking term is present into the lagrangian. In the 
following we study both cases. 

1. Chiral restoration in the chiral limit (e[ — 0) 

The restoration of chiral symmetry at hnite density 
and temperature, occurs as a first order transition only 
when the chiral invariance of the lagrangian is exact, or 
equivalently, only when the explicit symmetry breaking 
term = 0. In this case, at a given temperature we 
observe that for densities larger then a critical one the 
mean field value of a drops to zero and the masses of 
the chiral fields become degenerate. In Fig. [T] the behav- 
ior of the o~ field, as a function of the density and for 
various temperatures is displayed. The transition turns 
out to be first order, due to the discontinuous behavior 
of the sigma field. Similarly, in Figs. [H [3] we show the 
effective masses of the chiral fields, which become degen- 
erate when the symmetry is restored. The chiral phase 
diagram for isospin symmetric nuclear matter is shown 
in Fig.|¥]in the temperature versus baryon density plane. 
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FIG. 1: Mean value of the sigma field as a function of the baryon 
density for different values of the temperature. Here e' x = 0. 
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FIG. 2: Effective sigma meson mass as a function of the baryon 
density for different values of the temperature. Here e' x = 0. 



At zero baryon density the temperature at which the chi- 
ral restoration occurs is Tc = 165 MeV, which is rather 
close to the QCD critical temperature estimated by lat- 
tice QCD [HI H3- Obviously, in the case of QCD the 
restoration of chiral symmetry comes together with the 
appearance of the quark-gluon plasma phase. As it can 
be observed, the critical temperature decreases with in- 
creasing density, but it never reaches zero, not even at 
very large densities. Numerically it tends asymptotically 
to ~ 40 MeV. It can be interesting to notice that some 
(very preliminary) lattice calculations at finite density 
seem to support this result [28| . The chiral phase restora- 
tion in the chiral limit is also represented in Fig. [7j in the 
energy-density vs. baryon-density plane, by two closely 
spaced solid lines. Since the chiral transition is first order 
these lines are separated by a small energy gap. Indeed, 
at finite temperature a first order transition occurs with 
a jump in the energy density and in the entropy. 
The temperatures investigated in this Subsection are not 
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FIG. 3: Effective pion mass as a function of the baryon density 
for different values of the temperature. Here t' ± = 0. 
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FIG. 4: Chiral phase diagram for symmetric nuclear matter in 
the temperature versus baryon density plane. Here e\ = 0. 



to T ~ 196 MeV reaching very small but not vanishing 
values. This behavior, called "cross over", is predicted 
by QCD lattice calculations with two flavors at zero or 
small densities. 

At T > 196 MeV the dilaton field drops to zero (initially 
in the small density region, see Fig. [S]) and scale invari- 
ance (and therefore also chiral symmetry) is restored, as 
discussed in the next section. 

The large density behavior of the a mean field can be 
obtained by solving its equation of motion, what can be 
done analytically for T = 0. For large values of the scalar 
density ps, we obtain a oc 1/ps- Since ps scales as p 1 ^ 3 
for large p, the a mean field decreases monotonously but 
slowly with p and it never vanishes. This analytic result 
remains numerically true also at finite temperatures (for 
T < 196 MeV). 

Let us now concentrate on Fig. where a sort of 
"phase map" for symmetric nuclear matter is shown 
in the energy-density versus baryon-density (p, e) plane. 
The region in which the variation of the a field is more 
rapid is indicated by a dashed curve, which is obtained 
by joining the points of the (p, e) plane where the mod- 
ulus of the derivative of a, with respect to the temper- 
ature, takes its maximum value. In Ref. [29] the maxi- 
mum energy and baryon densities reached at the moment 
of maximum compression during a HIC experiment are 
computed within several models. It is interesting to no- 
tice that in all model estimates when the beam energy 
exceeds ~ 5 A GeV the maximum energy density lies 
above the dashed line of Fig. [7] and therefore signatures 
of the partial restoration of chiral symmetry should be 
observable. In particular, the softening of the EOS due 
to partial chiral symmetry restoration will be discussed 
later in this Section. 



large enough for the dilaton field to change significantly 
from its zero temperature value, therefore we could de- 
couple the discussion of chiral symmetry restoration from 
the restoration of scale invariance. This simplified sce- 
nario will not hold in the more realistic e[ =/= case, 
discussed in the next Subsection. 



2. Chiral restoration with a chirally broken 
lagrangian, e[ 7^ 

In the more realistic case in which t\ 7^ and the pion 
has a finite mass, we observe that the a mean field is 
continuous everywhere and it never drops to zero, at any 
density as long as T < 196 MeV, see Fig. [5] Therefore 
a chiral phase transition, as that observed in the chiral 
limit case, does not occur. On the other hand, the a mean 
field decreases continuously with the density and drops 
rapidly when the temperature raises from T ~ 100 MeV 
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symmetry is explicitly broken (e^ 7^ 0). 
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FIG. 6: Mean value of the dilaton field as a function of the baryon 
density for different values of the temperature. Here 7^ 0. 
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gap is of order of 1 GeV/fm 3 , a value not too different 
from the jump in energy density seen in lattice QCD in 
the pure gauge sector. Clearly, a direct comparison is 
meaningless since in full QCD the dcconfincment transi- 
tion (at zero baryon density) is a crossover and, moreover, 
we do not have quarks in our model. 

At zero baryon density, we obtain a critical temper- 
ature T c (p = 0) ~ 196 MeV (see Fig. |5J}, which corre- 
sponds to a jump in energy density from e~(p = 0) = 725 
MeV/fm 3 to e+(p = 0) = 1565 MeV/fm 3 . This critical 
temperature is significantly lower than the one estimated 
in H, 0] . This reduction is due to the contribution of the 
vector mesons, which provide a huge contribution to the 
pressure when the scale invariance is restored (since their 
masses drop to zero). For instance, in Ref. [|[, where no 
vector meson is present .the scale restoration occurs at 
T ~ 550 MeV. In Ref. [H the presence of the u> meson 
lowers the critical temperature to T ~ 300 MeV, while 
the inclusion in our calculations of the p meson (together 
with lu) lowers the restoration temperature to T ~ 200 
MeV. 

We notice from Fig. [7] that the critical energy density 
e c (p) increases with p. Also the critical temperature for 
scale invariance restoration increases with p, as shown in 
Fig. [5) This interesting effect arises from the repulsive 
contribution of the to vector meson, which contributes to 
the energy density as E u — (1/2) (^/m* ) 2 p 2 (here for 
simplicity we have assumed G4 = 0). Since the uj me- 
son mass scales with the dilaton field (see Eq. ([T5])), E u 
increases when the dilaton field decreases. Therefore at 
large p it is less convenient to restore scale symmetry [66[ . 
This effect is also clearly visible in Fig.|5]which shows how 
the dilaton field increases again at large densities. 
A few caveats are in order: 



FIG. 7: Phase diagram for symmetric nuclear matter in the 
energy-density versus baryon-density plane. Two closely spaced 
solid lines separate the (low-lying) region where chiral symmetry 
is broken from the (upper-lying) region where chiral symmetry is 
restored in the chiral limit = 0. All other graphical signs refer 
to the broken symmetry case t' x 7^ 0. The lower shaded area indi- 
cates the forbidden region under the T=0 EOS. The upper shaded 
area separates the (low-lying) region where scale symmetry is bro- 
ken from the region where scale symmetry is restored. The dashed 
line indicates the maximum (negative) variation of the cr-field mean 
value with the temperature (see text). 



3. Scale symmetry restoration 

In the model here discussed, the restoration of scale in- 
variance takes place when the mean value of the dilaton 
field drops to zero. As shown in Fig. [f)]the restoration of 
scale invariance takes place (in the mean field approxima- 
tion) as a first order transition, signaled by the discon- 
tinuous behavior of the dilaton field. A important result, 
shown in Fig. is that a large energy gap appears be- 
tween the two phases (due to the large amount of energy 
released by the condensates which drop to zero). This 



• we are working in a mean field approximation and 
it is well known that the order of a phase transition 
cannot be determined in that approximation; 
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FIG. 8: Phase diagram for symmetric nuclear matter in the tem- 
perature versus baryon density plane. 
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FIG. 9: Bulk modulus B = pdP/dp as a function of the baryon 
density at T=50 MeV for both the CDM and the Walecka GM3 
models. Here various values of the proton ratio are considered. 



• at the moment we have not included the thermal 
fluctuations of the dilaton field, while we have taken 
into account the thermal fluctuations of the vector 
mesons. The reason for this choice is that the con- 
tribution of the vector mesons largely dominates 
the contribution of the scalar dilaton field, since 
the contribution of the vector (isovector) fields is 
multiplied by their degeneracy. Moreover, to take 
into account the dilaton thermal contribution one 
has to overcome the difficulties associated with the 
logarithmic term. An expansion and resummation 
similar to the one performed for the chiral fields 
should be done, what we plan to do in a future 
work. 



B. Bulk modulus 

Here we discuss bulk properties of matter within the 
CDM. In particular we show in Fig. [5] the value of the 
bulk modulus B = pdP/dp as a function of p at T=50 
MeV. Here the calculations are performed for various 
isospin densities and we compare with the same quantity 
computed within the Walecka GM3 model [3(J. Notice 
that the bulk modulus computed with the CDM is larger 
than the one calculated with GM3 when p < 3 po. This is 
due to the larger value of the CDM incompressibility pa- 
rameter at po (K^ 1 = 9p 2 9 qJ± A — 322 MeV) compared 
with that obtained using GM3 (Ar 1 = 240 MeV). 

Observe that for p > 3po the effect of the partial 
restoration of the chiral symmetry manifests as a soften- 
ing of the EOS. This softening is clearly visible in Fig. [5] 
where the bulk modulus computed with the CDM equals 
the one evaluated using GM3 at about 3po and becomes 
smaller for larger values of p. The effect of the isospin 
density is to increase the bulk modulus, since the p meson 
gives a positive contribution to the pressure proportional 



to the baryon density. The increase of the bulk modulus 
when moving from Z/A=0.5 to Z/A=0.2 is of about 10% 
both for the CDM and for GM3. 

The dependence of the bulk modulus on the density 
and the temperature can be measured in HIC experi- 
ments, in particular by analyzing the so called collective 
flows. A first attempt has been done in Ref. [3l|, where 
they have extracted the value of the incompressibility 
parameter from HIC experiments, observing a significant 
softening of the EOS when the beam energy varies from 
2 A GeV up to 10 A GeV. By comparing with the CDM 
one can notice that this softening can be interpreted as 
due to partial chiral symmetry restoration if densities of 
at least 6 — 7po are reached during a HIC. In Ref. [13] we 
have shown that if such large densities are not reached, 
then this softening can better be interpreted as due to 
the formation of a mixed hadron-quark phase. 



C. Effective meson masses 

A very important point concerns the predictions of the 
CDM about the experimentally measurable "in medium" 
masses, as the vector meson masses and the pion mass, 
whose expressions are provided in Eqs. (|15I16I18[) . 

1. Effective pion mass 

The binding energy of the deeply bound Is and 2p 
states in pionic atoms of 207 Pb, established experimen- 
tally through the 207 Pb(d, 3 He) reaction, has been used 
to derive the value of the pion mass at nuclear matter 
density [32, 33] . The result of this still preliminary anal- 
ysis is that the pion mass increases by about 28 MeV at 
Po in symmetric nuclear matter. 

In Fig. 110b , the pion effective mass in symmetric nu- 
clear matter is plotted as a function of the baryon den- 
sity and for different values of the temperature. We ob- 
serve that the pion mass increases monotonously with 
the baryon density reaching a value of about 167 MeV 
at po and T = 0,a result very close to the one sug- 
gested in Refs.j32l.l33l]. We also compare the technique of 
Refs. @, |3] with the technique of Ref. [24[ (see the pre- 
vious discussion in subsection III C[) . Most of the results 
obtained in the present work are almost independent of 
which technique is adopted, but the masses of the chiral 
fields are very delicate and the two techniques provide 
different results even in the case of isospin symmetric 
matter, as we show in Fig. [TUb 

In Fig. [TUb we show the effect of the isospin asymme- 
try (Z/A=0.3) on the pion mass. Notice that in isospin 
asymmetric matter the positive and negative pions have 
opposite (nonzero) chemical potentials. Then, accord- 
ing to Eq. ([9]) , their thermal fluctuations are different as 
well as their effective masses, defined by Eq. (|T5)) . In 
order to consider this effect, in our calculations we used 
the technique of Ref. pij ]. In Fig. [TUb we show that the 
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FIG. 10: In panel (a) the pion effective mass obtained using the 
technique of Ref. [24J (MME, thick lines) is compared with the 
mass obtained using Refs. [3, 4] (CER, thin lines). The mass is 
plotted as a function of the baryon density, for Z/A=0.5 and for 
three values of the temperature. The two techniques are equivalent 
at T = 0. In panel (b) we show the masses of the different isospin 
components of the pion, as a function of the baryon density, for 
Z/A=0.3 and for two values of the temperature. Here we adopted 
the technique of Ref. |24| . At T = no isospin mass splitting is 
present. 



mass splitting of the different pion components increases 
substantially with the density and with the temperature 
(the mass splitting vanishes at T=0). In particular, at 
p = 3po and T=130 MeV, the difference between inl- 
and m OT + reaches a value of about 100 MeV. Notice that 
this result, obtained without including any contribution 
associated with the so-called particle-hole excitations, is 
comparable to the isospin splitting of the pion mass as- 
sociated with the so-called Tomozawa- Weinberg term. 



2. In medium vector meson masses 

The experimental determination of the in-medium 
modification of the light vector mesons (u>, p and $) 




FIG. 11: Vector meson masses as a function of the baryon density 
at T=0 and for different proton fractions. 



is still uncertain. Recent experiments have advanced 
to a level that in-medium spectral information on vec- 
tor mesons can be deduced. Experimentally, the best 
approach is the dilepton spectroscopy, in particular the 
investigation of the dileptonic decay modes of the light 
vector mesons, produced off nuclei or in HICs. In Fig.fTTI 
the masses of the vector mesons u> and p, evaluated within 
the CDM, are plotted as functions of the baryon density 
at zero temperature, and for different values of the proton 
ratio. We do not show the values of the vector mesons 
masses for temperatures different from zero because, up 
to T ~ 150 MeV, the masses vary only slightly with the 
temperature (the differences are of a few per cent). 

Notice the peculiar behavior of and m p with p: for 
low densities the masses initially drop, reaching a mini- 
mum in the neighbor of p — 3.5po, while for higher values 
of p, the masses increase again. This effect is due to the 
difficulty in restoring the scale invariance at large densi- 
ties, as discussed in subsection MI A 31 The special sce- 
nario obtained within the CDM seems to be in agreement 
with several experimental analysis indicating a moderate 
reduction of the masses of the vector mesons at densities 
of the order of po and no significant reduction when the 
masses are tested at very large densities. 

For instance, in experiments of scattering of photons 
and protons on nuclei, densities about po are tested and 
the results indicate a drop of the w-mass by 14% at satu- 
ration density (CBELSA/TAPS) [II]. More controversial 
is the situation concerning the p-mass: the KEK-PS E325 
collaboration indicates a drop by 9% at p [35|, HH while 
a recent experiment at JLAB suggests a smaller drop by 
~ 2% at po [37|. On the other hand, in HIC experiments 
at energies of the order of one hundred A GeV (SPS and 
NA60), the vector mesons mass shift has been found to 
vanish [36|, [H, H§| . This result presumably discards the 
scenario predicted by Brown and Rho [36l . |40| | , where the 
vector meson masses drop continuously with the density 
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Finally, the effect of a finite isospin density is to in- 
crease the vector masses as it can be observed in Fig. [TT] 
This effect grows with the baryon density, although at 
p = 9po the increment is only of about 1% for both the 
vector meson masses. 



IV. ASTROPHYSICAL APPLICATIONS 

Neutron Stars (NSs) are the most dense objects in na- 
ture while their temperature is always below a few ten 
MeV and typically lower than one MeV. Therefore they 
are the only objects which can provide information on 
the EOS of matter at very large densities and low tem- 
perature. Unfortunately the structure of compact stars is 
still uncertain since astrophysical data can provide only 
rather indirect indications, leaving the field open to theo- 
retical investigations based on models for the high density 
EOS. Since the CDM is a chiral model, it is very interest- 
ing to investigate its previsions concerning matter at very 
large densities and small temperatures. In particular it 
is possible to explore how the chiral dynamics can affect 
the structure and the formation of a NS. The dynamics 
of the gluon condensate does not play a crucial role at 
low temperatures. In this section we compute the value 
of the adiabatic index which plays a crucial role in Su- 
pernova explosions and we investigate within the CDM 
the structure of NSs [lj| • 

In order to better understand the formation and the 
structure of a NS we also discuss the effective masses of 
the various fields and, finally, we consider the possibility 
to reach the pion condensation in NSs. 



A. Adiabatic index 

As far as we know, the explosion of a core-collapse su- 
pernova proceeds as the following (for a recent review 
see Ref. (4l|): when a massive star exhausts its nuclear 
fuel, the gravitational force cannot be contrasted by the 
outward pressure and the core starts to collapse. During 
this stage, which typically last much less than a second, 
the pressure strongly exceeds the electron degeneration 
pressure and neutronization occurs. The temperature is 
of the order of few ten MeV and the neutrinos cannot es- 
cape from the collapsing matter since their cross section 
is very large. The lepton number is then conserved and 
no energy can escape the system, giving the possibility to 
describe the entire process as an adiabatic compression 
with a entropy per baryon of the order of unity. More- 
over, since the time scale of the (3 interactions is of the 
order of 1CP 8 s, the system has the time to reach the 
/3-equilibrium. When the density of the core overcomes 
po, due to the nuclear repulsion the core bounces trig- 
gering a shock wave which propagates to the outer layers 
of the star. In hydrodynamical simulations, it turns out 
that the shock wave stalls while it is still propagating 



inside the outer layers of the progenitor star. Therefore 
the so called "direct mechanism", characterized by the 
idea that the explosion is directly related to the shock 
wave generated by the bounce, has been ruled out. The 
failure of this direct mechanism sets the stage for a 'de- 
layed' mechanism, in which the shock is re-energized by 
the heating of the outer layers due to the intense neutrino 
flux emerging from the neutrino-sphere. 

It is worth discussing more in details the conditions 
under which the direct mechanism would be successful. 
It is has been noticed that the possibility for a supernova 
to explode via the "direct mechanism" is related to the 
softening of the EOS at densities just above po (42J. This 
possibility, exploited in the Baron-Cooperstein-Kahana 
(BCK) EOS [43| has been discarded because it seems to 
be inconsistent with the constraints on the NS masses 
coming from observations, since a too soft EOS cannot 
support a sufficiently massive NS. A more sophisticated 
approach invokes a softening taking place only in a finite 
density window and associated with a phase transition. 
In this way the EOS can be stiff again at densities larger 
than the ones associated with the transition. This possi- 
bility is also not new and it has been investigated in sev- 
eral papers [3, [45|, |H, [4?J ■ In all those works a strong 
softening of the EOS is ascribed to a phase transition 
which can be due either to the pion condensation or to 
the formation of a mixed phase of hadrons and quarks. 
Since in the CDM we obtain a softening of the EOS at 
moderate densities due to chiral symmetry restoration it 
is worth to explore the possible influence of this transi- 
tion on the supernova explosion. 

In order to investigate the core collapse of a massive 
star, we compute the EOS along adiabatic paths for mat- 
ter composed of neutrons, protons, electrons and elec- 
tronic neutrinos. Due to the previous considerations, the 
conditions imposed to the system are: 

i) leptonic number conservation: 

Pe + p Ve = Y ia (p p + p n ) (22) 

ii) electric charge neutrality: 

P P = Pe (23) 

iii) /3-equilibrium: 

Pn - Pp = Pe- Pv, (24) 

where Yj c = 0.4 is the constant ratio of the leptonic den- 
sity to the baryon density. Since before the collapse no 
muons are present, we keep — Pll+ ^ = 0. 

An important quantity which allows to estimate the 
softness of the EOS is the adiabatic index, which is de- 
fined as: 

r = §^ (25) 

amp 
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FIG. 12: Adiabatic index as a function of p computed within 
various models for isoentropic n,p, e~,v e matter with Yj = 0.4. 
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FIG. 13: Pressure for NS matter as a function of baryon density, 
evaluated within the CDM and within the Walecka GM3 model. 



where P is the pressure. It corresponds to the index 

of a polytropic equation of state P(p) = Pa (~^J • ln 
Fig. [12] we show T, obtained from the CDM EOS and 
from other models, as a function of p. Notice that within 
the Walecka GM3 and the Lattimer-Swesty EOS [H[ no 
significant reduction of T is observed at supranuclear den- 
sities. At the contrar y, t he strong softening occurring 
within the BCK EOS [43| at densities just above po al- 
lows a supernova to explode via the "direct mechanism" 
(recall that the BCK EOS has been ruled out since it 
provides too low NS masses). Within the CDM the par- 
tial restoration of chiral symmetry manifests as a strong 
softening of the EOS (see Fig. IT5|) . This softening is also 
observed as a strong reduction of T for densities larger 
than 3po- Although this reduction is comparable with 
that obtained in BCK, since the density at which it oc- 
curs is much larger than po we do not expect that the 
CDM can allow a supernova to explode via the "direct 
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FIG. 14: (Color online) Mass-radius relation for the CDM and for 
various other models together with recent experimental constraints. 
Adapted from Ref. 



mechanism" . On the other hand, since very large densi- 
ties are reached during the collapse of a binary system of 
neutron stars, we expect that such a strong decrease of 
r can have interesting consequences, for instance in af- 
fecting the gravitational wave signal emitted just before 
the formation of the black- hole I50| . 



B. Neutron stars 

A NS is the cold remnant of the explosion of a core- 
collapse supernova and it is composed by neutrons, pro- 
tons, electrons and muons in /3-equilibrium. Since the 
temperature is very low the neutrinos cross section is 
very small, allowing the neutrinos produced by /3-decay 
to escape from the star. Therefore the lepton number is 
not conserved. Since electric charge neutrality has also 
to be imposed, the conditions defining NS matter read: 



Pn 
Pp 



= p p + p e 
Pe+ Pp- 



(26) 



That EOS is then used to solve the Tolman- 
Oppenheimer-Volkov (TOV) equation obtaining a se- 
quence of stable stellar configuration. 

The mass-radius relations for the CDM and various 
other models are shown in Fig. 1 141 together with sev- 
eral recent experimental constraints coming from obser- 
vations and discussed in Ref. [49( . The models considered 
in Fig. Q3] can be separated into Relativistic Mean Field 
(RMF) phenomenological models, including the Walecka 
type models (GM3 |3, NLp [5l| and NLp5 [H) and 
those with density dependent coupling constants (DD 
[H, D 3 C andDD-F ^Si), the "ab-initio" calculations 
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FIG. 15: (Color online) NS mass vs. central density for the same 
models of Fig. 1141 The full circles indicate the DU thresholds and 
the stars indicate the maximum mass configurations. Adapted from 
Ref. 



including the relativistic Dirac-Brueckner-Hartree-Fock 
(DBHF MpLilllp) and the hybrid models (KVR and 
KVOR [5y, |60j) which bridge the gap between the fully 
microscopic and more phenomenological descriptions. 

The experimental constraints are extracted from obser- 
vations of kilohcrtz quasi-periodic brightness oscillations 
in the Low Mass X-Binaries 4U 0614+09 and 4U 1636- 
536 and from the thermal radiation of the isolated NS 
source RX J1856.5-3754. 

Due to the rather unrealistically large value of the in- 
compressibility parameter at po (K^ 1 —322 MeV), larger 
NS radii are obtained within the CDM than within the 
other models. Therefore, although the CDM allows to 
have a NS with a mass M ~ 1.4M Q and a radius of 
14 Km, satisfying the constraints from RX J1856 (the 
other models exclude this interpretation), our result is 
questionable due to the too large value of the incom- 
pressibility at saturation. 

An important check is based on the recent observations 
of PSR J0751+1807. The data analysis, although still un- 
certain, allows to extract a NS mass of 2.1±0.2 ftois) Me 
(first error estimate with la confidence, second in brack- 
ets with 2a confidence) [6l|. Since the CDM exhibits a 
maximum mass M — 2.06M© this constraint is easily 
fulfilled. The constraint is displayed in Fig. [T5l where the 
NS mass vs. central density relation is shown for the 
same models of Fig. [14] and where the crosses denote the 
maximum mass configurations. 

When discussing NSs a very important role is played by 
the cooling mechanism. We can roughly divide the cool- 
ing mechanisms in "slow" ones, as e.g. the so-called mod- 
ified Urea process, and fast ones, as e.g. the Direct Urea 
(DU) process which corresponds to a sequence of neu- 
tron /3-decays and lepton captures (see e.g. Ref. ([62|)). 
From simple considerations it is possible to find the crit- 



ical baryon density beyond which the DU processes can 
occur in a NS, which turns out to be the density at which 
the proton ratio in NS matter reaches the value: 



where x e = p ^ p . In Fig. [15] the DU-thresholds are de- 
noted by full circles. Notice that, although the CDM al- 
lows the DU to occur for NS masses exceeding ~ 1.2M©, 
its threshold is slightly higher than those obtained within 
GM3, NLp, NLpS and DBHF. 

Finally, the effect of the softening of the CDM EOS 
manifests as an "anomalous" behavior of the NS sequence 
with respect to the other models. Notice in fact that 
in Fig. [15] at low central densities the CDM NS masses 
are rather large and the CDM sequence seems to follow 
the trend of the DD, D 3 C and DBHF models whose 
maximum mass configurations exceed 2.3-M©. However 
at larger central densities the sequence flattens and a 
maximum mass only slightly larger than GM3 is reached 
(M max = 2.06M Q ). 

C. In medium masses 

In order to understand how the partial restoration of 
chiral symmetry affects the structure of NSs, we com- 
puted, within the CDM and using the technique devel- 
oped in Ref. (24|, the n° and a meson effective masses 
as a function of p both for NS matter and for isoentropic 
neutrino-trapped matter (see Fig. \W\i . We find that even 
at large densities the masses of the two chiral fields are 
still rather different, although the mass of ir° is consid- 
erably larger than in the vacuum. The effect of the tem- 
perature on the partial restoration of chiral symmetry is 
very small but it can anyway be noticed by comparing 
the isoentropic curves with S/R=l and S/R=2: for the 
larger value of the entropy the a and tt° masses start get- 
ting closer, the a mass by decreasing and the ir° mass by 
increasing its value. 

Finally, we have examined the in-medium modification 
of the nucleon and of the pion mass by plotting their 
values as a function of the radius in the case of a spherical 
non-rotating NS of mass M = 1.4M©, see Fig. [T71 (notice 
that, although in presence of isospin asymmetry, at T=0 
it comes out that m^o = m„+ = m^- = m-rc). The value 
of nucleon mass decreases by more than 50% when going 
from the outer crust to the center of the star. The pion 
mass, at the contrary, increases by about 45%. 

D. Pion condensate 

Since in compact stars very large baryon densities and 
isospin asymmetries are reached, it is interesting to in- 
vestigate the possibility to form a meson Bose-Einstein 
Condensate (BEC) . In particular we investigate here the 
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FIG. 16: a (figure (a)) and 7T° (figure (b)) effective masses as a 
function of the baryon density for NS matter (dashed line) and 
for isoentropic ^-trapped matter with entropy per nucleon S/R=l 
(thin solid line) and S/R=2 (bold solid line). 
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FIG. 17: In-medium mass modification of the nucleon (figure (a)) 
and of the pion (figure (b)) inside a neutron star of 1.4 Mq. 



possibility of reaching the pion BEC, which is easier to 
achieve due to the small value of the pion mass. Usually 
with "pion condensate" one refers to a collective pionic 
excitation at non zero-momentum. This collective state 
is a pole in the pion propagator modified by the inclu- 
sion of the particle-hole excitations. In this work we do 
not include the particle-hole excitations in order to sat- 
isfy the Goldstone theorem at finite temperature (see the 
discussion in section Hi E[) and we call "pion condensate" 
the zero-momentum Bose-Einstein condensate of pions. 

Recall that in order to conserve the isospin charge a 
chemical potential is assigned to the pions, which reads: 
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FIG. 18: 7r~ effective mass and effective chemical potential as 
a function of the baryon density, computed within the CDM for 
NS matter and for isoentropic ^-trapped matter with entropy per 
baryon S/R=l. 



where p p and p n are the chemical potential of protons 
and neutrons, respectively. The equation above assures 
the equilibrium in the reaction n <-> p + n~ . Moreover 
the pion effective chemical potential, which enters the 
thermodynamical integrals, is given by: 



(29) 



where p* and p* n are the effective chemical potential of 
protons and neutrons respectively, given by Eq. (|19p . The 
effective chemical potential of the pion increases with the 
baryon density and with the isospin asymmetry. When 
the effective chemical potential becomes equal to the pion 
effective mass, real pions start to be produced in the 
medium and they condense in the state of minimum en- 
ergy, i.e. the zero momentum mode. The condition for 
the onset of the pion BEC then reads [25j : 



(30) 



/i 7r + j-ip p n , 



(28) 



In Fig. [18] the ir~ effective mass and its effective chem- 
ical potential (computed using the technique developed 
in Ref. are plotted as a function of p. Here we 

consider two different cases: /3-stable NS matter at T=0 
and isoentropic z^-trapped matter at /3-equilibrium with 
entropy per baryon S/R=l. As it can be seen the ir~ 
mass grows with p and it does not change appreciably 
between the two cases. Also fi v - increases with p, but 
its value actually depends on the value of the entropy 
and on neutrino-trapping. 

The important result is that the condition for the onset 
of the BEC is never achieved, not even at large densities. 
This is due to the increase of the pion effective mass with 
p, which remains always larger than the chemical poten- 
tial. In simpler models, like e.g. the Walecka GM3 where 
the pion mass is constant, the increase of the chemical 
potential gives rise to the pion condensation, as discussed 
in Ref. 1631. 
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V. CONCLUSIONS 

The main aim of this work was to investigate the be- 
havior of matter at large densities and temperatures by 
using a effective lagrangian with broken scale invariance 
and chiral symmetry. Our results are concordant with 
those obtained in Refs. HQ. We extended the previous 
calculations in several directions and wherever possible 
we compared with experimental data. We have also used 
two different techniques to evaluate the effects of thermal 
fluctuations, i.e. the method developed in Refs. [3, B| 
versus the one developed in Ref. [24(. The latter is more 
general since it does not make any assumption concerning 
the value of the thermal fluctuations of the chiral fields. 
The two techniques provide similar results as long as we 
are not investigating subtle isospin dependencies such as 
the mass splitting of the isospin components of the pion. 

Our main results are the following: 

- we have provided a phase diagram mapping the regions 
in which chiral symmetry and/or scale invariance are re- 
stored. The scenario provided by the model we have 
studied, i.e. the decoupling of chiral symmetry restora- 
tion and scale invariance restoration (which is likely as- 
sociated with deconfinement) shares strong similarities 
with the scenario recently suggested by McLerran and 
Pisarski in the large N c limit [641 ]; 

- the masses of vector mesons are affected by the partial 
restoration of scale invariance, since they scale with the 
dilaton field. Therefore scale invariance restoration is 
more difficult at large densities, because the repulsive 
effect of the exchange of vector mesons is enhanced if 
their mass is reduced; 

- the masses of the vector mesons first reduce at finite 
density, but at larger densities they even increase. This 
result holds up to temperatures of the order of 200 MeV, 
at even larger temperatures scale invariance starts being 
restored and the vector meson masses drop to zero; 

In isospin symmetric matter the pion mass increases 



with the density and with the temperature. At po and 
T=0 the effective mass of the pion is about 170 MeV. 
In presence of isospin asymmetry the masses of the differ- 
ent pion components split. The splitting increases with 
the density and the temperature and, for instance, at 
p = 3p , T=130 MeV and Z/A=0.3 m w - - m T+ ~ 100 
MeV; 

- a very substantial reduction of the adiabatic index takes 
place at densities exceeding 3 po- Although this softening 
cannot directly affect the mechanism of Supernova explo- 
sion, because it starts at too large densities, its signature 
could be observed during a two neutron star merging; 

- the BEC of pions does not form in NSs due to the 
increase of the pion mass with the density. 

The present work can be extended in several directions. 
In particular the restoration of scale invariance will be 
better investigated by including in the calculation also 
the thermal contribution of the dilaton field, while in 
the present work we have only taken into account the 
mean field effect. Hyperons should also be included in 
the model. This direction is now explored by analysis 
similar to ours and performed by the Frankfurt group 
[ilil ]. The most interesting extension would be to take 
into account quark degrees of freedom. An attempt in 
that direction is a recent paper in which deconfinement 
is discussed in a dilaton model, but in the pure gauge 
sector [Hj]. 
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